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ABSTRACT 

S-  '.  Asymptotic  solutions  are  given  for  temperature  profiles  in  laminar  three- 

dimensional  flows  with  viscous  dissipation.  The  results  are  asymptotically 

f' 

valid  for  small  thermal  dif  fusivity  ■  fi#'  they  hold  for  Newtonian  or  non- 
Newtonian  fluids.  The  fluid  properties  are  evaluated  at  an  average  temperature 
for  the  system;  this  is  satisfactory  for  moderate  temperature  differences. 

Heat  transfer  formulas  for  various  thermal  boundary  conditions  are  included. 
The  general  formulas  are  applied  to  several  geometries;  numerical  results  are 
obtained  for  a  wire-coating  operati'-..,  '■'xy' 
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SIGNIFICANCE  AND  EXPLANATION 


In  this  paper  we  study  the  effect  of  viscous  heating  on  the  temperature 
distribution  in  a  flowing  liquid.  Matched  asymptotic  solutions  are  obtained 
for  the  temperature  distribution  in  the  presence  of  viscous  dissipation,  with 
and  without  interfacial  heat  transfer.  The  results  are  asymptotically  valid 
for  small  thermal  diffusivity  (large  Prandtl  number);  they  should  prove  useful 
in  studies  of  wire  coating,  extrusion,  heat  transfer,  and  other  flow  processes 
with  high  rates  of  shear  or  highly  viscous  liquids 


The  responsibility  for  the  wording  and  views  expressed  m  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 
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FORCED  CONVECTION  IN  THREE-DIMENSIONAL  FLOWS:  III. 

ASYMPTOTIC  SOLUTIONS  WITH  VISCOUS  HEATING 

Matthew  A.  McClelland  and  Warren  E.  Stewart 

SCOPE 

Small-diffusivity  asymptotes  have  given  useful  information  on  a  variety 
of  heat  and  mass  transfer  systems.  Part  I  and  II  of  this  series  (Stewart, 
1963;  Stewart,  Angelo  and  Lightfoot,  1970)  dealt  with  boundary-layer  phenomena 
in  the  absence  of  viscous  heating.  The  present  paper  adds  the  effects  of 
viscous  heating  in  the  boundary  layers  and  in  the  main  flow  region.  The 
results  should  prove  useful  in  studies  of  wire-coating,  extrusion,  heat 
transfer,  and  other  flow  processes  involving  rapid  deformation  or  highly 
viscous  fluids. 

This  treatment  holds  for  laminar  non-separated  regions  of  external  flows, 
and  for  laminar  thermal  entrance  regions  of  internal  flows.  The  velocity 
profiles  are  considered  as  given  or  separately  calculable;  this  is  consistent 
with  the  use  of  an  average  temperature  in  evaluating  the  fluid  properties. 


Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-0041 


CONCLUSIONS 


The  key  results  of  this  analysis  are  the  inner  and  outer  solutions,  given 
in  Equations  (46),  (15),  and  (54).  These  may  be  combined  with  the 
superposition  formulas  of  Papers  I  and  II  (Stewart,  1963;  Stewart,  Angelo  and 
Lightfoot,  1970)  to  predict  the  temperature  profiles  and  heat  transfer  rates 
in  various  three-dimensional  flows  once  the  velocity  profiles  are  given.  The 
resulting  general  formulas  for  steady-state  systems  are  given  in  Eqs.  (47)  - 
(51 ). 

For  fluids  with  small  thermal  dif fusivities ,  a,  the  enthalpy  increase 

caused  by  viscous  heating  is  of  order  a  in  boundary  layers  that  are 

steady  as  viewed  from  the  nearest  interfacial  element.  The  increase  is  of 

order  in  time-dependent  boundary  layers,  and  is  of  order  o°  outside 

-  1/3 

the  boundary  layers  both  for  steady  and  unsteady  flows.  The  order  a  ' 
agrees  with  that  found  by  Meksyn  (1960)  for  the  steady  laminar  boundary  layer 
on  a  flat  plate  at  high  Prandtl  numbers.  A  somewhat  stronger 

-Vo 

dependence,  a  * ,  holds  for  steady  laminar  boundary  layers  at  Prandtl 
numbers  near  unity  (Pohlhausen,  1921). 

INTRODUCTION 

For  complex  flows  with  viscous  dissipation,  it  is  difficult  and  time- 
consuming  to  solve  the  equations  of  change  in  full  detail.  However,  useful 
asymptotes  are  obtainable  for  heat  transfer  by  considering  the  limiting  case 
of  small  thermal  diffusivity. 

Paper  I  (Stewart,  1963)  gave  asymptotic  solutions  for  forced  convection 
near  fixed  interfaces,  with  viscous  dissipation  neglected.  The  boundary 
layers  considered  there  were  at  steady  state  when  viewed  from  any  interfacial 


element 


Paper  II  (Stewart  et  al.,  1970)  dealt  with  forced  convection  near  mobile 
interfaces.  The  boundary  layers  considered  there  were  time-dependent  as 
viewed  from  the  nearest  interfacial  element.  The  analysis  included  internal 
heat  sources,  but  viscous  dissipation  was  not  explicitly  treated. 

In  this  work  the  results  of  Papers  I  and  II  are  extended  to  include 
viscous  dissipation.  We  use  the  continuity  and  energy  equations  for  a 
constant-property  fluid 

(V  •  v)  =  0  (1) 


DT  n2 
—  -aVT=- 


(£: Vv) 


PC 


(2) 


and  treat  the  velocity  field  as  given. 

Particular  integrals  of  Equation  (2)  are  derived  here  as  matched 
asymptotic  solutions  for  the  boundary  layers  and  the  main  flow  region.  The 
complementary  solutions  from  Papers  I  and  II  can  be  superimposed  on  these 
results,  to  describe  the  combined  effects  of  interfacial  heat  transfer  and 
viscous  dissipation.  We  do  the  superposition  explicitly  for  steady-state 


systems  in  Eqs.  (47)-(51) 


COORDINATES  AND  VELOCITIES 


The  steady-state  inner  (boundary  layer)  solution  is  done  in  the 
curvilinear  coordinates  x,  y,  and  z  of  Paper  I.  The  coordinate  x  is 
measured  along  the  surface  streamlines;  simplifications  for  two-dimensional 
systems  will  be  considered  later.  The  coordinate  y  is  the  normal  distance 
from  the  surface,  and  the  surface  coordinate  z  is  measured  normal  to  the 
local  x-direction.  In  the  neighborhood  of  a  given  interface  these  coordinates 
are  orthogonal,  with  scale  factors  defined  by 

dr=6hdx+6hdy+6hdz.  (3) 

~  ~x  x  ~yy  ~zz 

The  scale  factor  hy  is  unity  here,  and  hx  and  hz  are  taken  at  their 
inter facial  values. 

The  steady-state  outer  (main  flow)  solution  is  done  in  streamline 
coordinates,  illustrated  in  Figure  1.  Each  streamline  is  given  by  the 
intersection  of  two  stream  sheets  of  constant  ^  and  ^  ,  and  each  point 
on  a  streamline  is  identified  by  i  ,  the  arc  length  downstream  from  a 
reference  surface. 

The  velocity  components  for  the  steady-state  inner  solution  are  taken 
from  Paper  I.  They  reduce  to 

vxl}  -  y  B(x,z)  +  0(y2)  (4) 


ihiri7  lh.“>  * 0(y3> 

X  z 


(5) 


v(l)  =0  +  0(y2)  (6) 

z 

in  the  absence  of  net  mass  transfer.  These  expressions  satisfy  Equation  (1) 
to  first  order  in  y  .  The  non-negative  function  3  is  the  interfacial 
magnitude  of  the  rate -of -strain  tensor  as  defined  by  Bird,  Armstrong  and 
Hassager  (1977);  for  brevity  we  call  this  quantity  the  interfacial  shear 


Coordinates  for  outer  solution. 


rate.  A  particular  flow  field  is  described  (to  the  indicated  order  in  y)  by 
specifyinq  the  surface  coordinate  qrid  and  the  function  B(x,z). 


For  the  subsequent  matchinq  of  inner  and  outer  solutions,  we  need  to 
express  Equation  (4)  in  the  coordinates  of  Fiqure  1.  The  continuity  equation 
is  automatically  satisfied  by  settinq 

V  =  (V  ^  *  V  (7) 

which  may  be  applied  in  the  boundary  layer  coordinates  (x,y,z)  to  qive 


5 

~x 


V 


(i) 


h  3x 


This  expression  aqrees  with  Equations 
functions  in  the  forms 


6 

~y 


h  3y 


h  3y 


(4),  (5), 


5 

~z 


h  3z 
z 


(8) 


and  (6)  if  we  take  the  stream 


*  =  h  8  +  0(y3) 

I  z  Z 

♦2  =  z  +  C  +  0(y) 

with  C  an  arbitrary  constant.  Hence,  Equation  (4)  can  be  rewritten 

v(i)=  /TJ~T/h~  +  0(*_)  . 

x  1  z  1 


(9) 

(10) 


(11) 
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OUTER  SOLUTION  FOR  STEADY  FLOWS 


At  a  distance  from  the  boundaries,  the  heat  conduction  term  of  Equation 
(2)  can  be  neqlected  as  a  *  0  .  For  steady-state  conditions  this  equation 
then  becomes 


PC  (V  •  VT*°* )  =  -(t  :  Vv) 

p  **+  *w  <v 


which  may  be  rewritten  in  the  streamline  coordinates  as 

Jo) 

*  — (T  :  Vv)  . 


3t 

PCP  vsr 


Integration  of  this  equation  with  the  initial  condition 

T  =  T  at  1  =  0 

GO 

yields  the  outer  temperature  profile 

.  .  ,  ,  -(x  s  Vv) 

Tl  '(*,*»*,)  -  T.  -  -i-  /  - ^ -  d* 

1  1  PC  0  v 

P 


1 


(12) 


(13) 


(14) 


(15) 


which  is  well-behaved,  except  near  stationary  boundaries  where  V  vanishes 

and  (T;Vy)  does  not.  Near  such  boundaries,  heat  conduction  becomes  important 

and  an  inner  solution  is  required. 

The  subsequent  matching  of  solutions  will  require  the  asymptote  of 

Equation  (15)  as  ^  +  0  .  To  order  /i|’1  we  may  set  V  =  V^  * , 

dt  =  h  dx  ,  and 
x 

-  (r  :  Vy)  =  T  (x,z)8  (16) 

~  o 

in  which  tq(x,z)  is  the  wall  shear  stress.  The  resulting  asymptote  of 
Equation  (14)  near  stationary  boundaries  is 

lim  SW.  [T(o)  (*,♦  ,♦  )  -  TJ  =  -4r-  JX  T  /iTe  h  dx  (17) 

♦  ,♦0  1  12  PC  0  °  z  x  1 

1  p 

with  l  =  fx  h  dx.  in  this  limit. 

0  x  ’ 
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INNER  SOLUTION  FOR  STEADY  FLOWS 

Equation  (2)  can  be  simplified  according  to  the  usual  boundary  layer 

approximations  (Schlichtinq,  1960;  Rosenhead,  1963).  In  the  limit  of  small 

2 

o  ,  the  Laplacian  V  t  reduces  essentially  to  the  second  normal  derivative, 
82T 

— —  .  In  the  same  limit.  Equation  (16)  describes  the  viscous  dissipation  in 

3y 

the  boundary  layer.  With  these  substitutions,  and  use  of  Equations  (4),  (5) 
and  (6), the  enerqy  equation  becomes 


y£  ^T 

h 


(i) 


v  3 

—r-  a-  <h  0) 
2h  h  3x  z 
X  z 


3t 


3y 


(i)  32T(i)  t0p 

=  «  — ~  + 

3y^  PC 


(18) 


for  steady-state  boundary  layers  in  the  limit  of  small  a.  Since  no  z- 
derivatives  remain,  Ecniation  (18)  can  be  integrated  in  two  dimensions  with 
z  as  a  parameter. 

Eq.  (18)  can  be  solved  conveniently  by  representing  x^B  as  a 
superposition  of  special  functions  e(x,£,z).  Each  e-function  has  a 
different  value  of  the  location  parameter  £,  and  vanishes  for  x  <  £.  A 
corresponding  family  of  solutions  T*^  is  defined  by  the  differential 
equation 


yB  3t 


*(i) 


h  3x 
x 


y2  3  3T*(i) 

---  a~  (h_B)  1 - 


2h  h  3x  'V'  3y 
X  z 


32T*(1)  e(x,?,z) 

=  O  - -  +  - x - 


3y 


PC 


(19) 


and  the  boundary  conditions 


3T 


(i) 


3  v 


at  y  = 


(20) 


*  ( i  )  „ 

T  T_ 


as 


(21) 


*(i> 

T  =  T_ 


for  x  <  5 


(22) 
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Once  the  functions  e(x,£,z)  are  found,  these  and  the  resultinq  solutions 
T*^  can  he  superimposed  to  solve  Equation  (18)  for  an  adiabatic  boundary. 

*  t  i  ) 

We  postulate  the  followinq  form  for  T  , 


*  ( i  ) 

T  '  ’  -  T  =  f (x,£,z)q(n) 


(23) 


(24) 


with 

n  “  5(x,£,z)  * 

Here  S(x,S,z)  is  a  boundary  layer  thickness,  which  [in  view  of  Eq.  (22)]  is 
taken  to  vanish  for  x  <  £•  The  boundary  conditions  (20)  and  (21)  then  reduce 
to 


dq 


dn 


=  0  at  n  =  0 


(25) 

g  +  0  as  h  ♦  “  .  (26) 

Insertion  of  the  postulates  (23)  and  (24)  into  Equation  (19)  gives 


&£]*-[  ii  £ . 

X 


T  «L_  «  (h  8)]n2  £  -  *^2  +  [  ]  .  (27) 

ah  3x  2ah  h  3x  z  J  dn  ,_2  kf  J 


x  z 


dn 


The  postulates  hold  only  if  the  bracketted  terms  are  independent  of  x  ,  for 
x  >  (  ; 

£  £  ♦  i-  (h  6, 

ah  3x  2ah  h  3x  z 
X  X  z 


i!/n  , 

«nr-  17  <s  /s7>  -  <V*> 

X  z 


(28) 


6  j5  3  Inf 

T 

X 


5h“  ~5T~  =  ci(z) 


ST"  C2(Z)* 


(29) 


(30) 


The  functions  C^(z)  may  be  chosen  as  positive  constants.  Choosinq  =  3 


and  inteqratinq  Equation  (28)  with  the  initial  condition 


'h  3  S(x,C,z)  =  0 

7, 


at  x  =  S 


(31  ) 


required  by  Eq.  (22),  we  obtain  the  boundary  layer  thickness 

.1/3 


5  =  — - —  {9a  /X  /h  8  h  h  dx  } 
r. — s-  i  z  x  z  1 


(x  >  £). 


(32) 


'h  6 
z 


The  choice  C2  =  1  qives 

f  =  62c/k  .  (33) 

Equations  (28)  and  (29)  qive  the  followinq  differential  equation  for  f  : 


3*n  f 
3x 


C  3 Jin  (  5/h  $) 
1  _ z _ 

c  5x 

0 


(x  >  £ ) 


(34) 


Let  c-|/c0  =  r  =  constant;  then  Equation  (34)  may  be  integrated  to  qive 

f  =  C  (5  /hl)r  (x  >  £).  (35) 

3  z  — 

The  value  of  is  determined  by  matching  the  inner  and  outer  solutions 

for  the  reqion  x  >  £.  The  outer  limit  of  the  inner  solution  has  the  form 

U.  (HT‘lil  -  T.l)  ■  '  a  06) 

n-*»  0  i 

according  to  Equations  (23)  and  (27)-(30).  Insertion  of  Equation  (35)  gives 


lim  {nrrMl>  -  T,,,]} 
n-M» 


C  5r(h  0)r/2 


(37) 


The  inner  limit.  Equation  (17),  of  the  outer  solution  may  be  rewritten  for  the 
source  distribution  e(x,£,z)  to  give 


lim  KfW  (T*(o)  -  Tm)}  =  -i-  /X  e/jTTP  hxdx1  •  (3B) 

^♦0  1  PCp  £ 

Here  £  has  replaced  zero  as  the  lower  limit  of  integration.  Insertion  of 

Equations  (9),  (33)  and  (35)  qives  r_, 

C  01  1 

lim  {yfT*(o)-  T  ) )  =  r  h  h  «r'2(h  B)  2  dx  .  (39) 

y+0  /O  £  x  2 

z 
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For  r 


2 


and  C1  =  6,  Equations  (37)  and  (39)  match, 

lim  {y[T*(l)  -  T  J }  =  lim  (y[T*(o)  -  T—] } 
y+«  y+0 


giving 


cisV 


(40) 


Equation  (27)  then  becomes 

eng  -  3h2  -  ^-2-  +1  (41) 

dn  dn2 

which  may  be  solved  numerically  with  boundary  conditions  (25)  and  (26).  The 
resulting  surface  value  of  the  function  g  is 

g(0)  =  0.32973  (42) 

and  the  function  g(n)  is  shown  in  Figure  2. 

The  constant  of  Equation  (35)  remains  to  be  determined.  Equations 

(33)  and  (35)  give,  for  r  =  2, 

C3  =  kK <x>£).  (43) 

z 

By  definition,  e  vanishes  for  x  <  £;  thus,  e/h  0  is  a  step  function. 

z 

Insertion  of  Equations  (43),  (32),  and  (35)  with  r  =  2  into  Equation  (23) 
then  gives 


T*(l)  -  T  = 


>2/3 


9(h)  l  /X  -fiTS  hxhzdXlJ2/3  —  (44) 

-  e  z 


PC  .V3  x  =£ 

P  1 


as  the  fundamental  inner  solution. 

We  now  return  to  the  solution  of  Equation  (18)  with  an  adiabatic  boundary 
at  y  =  0  ,  a  uniform  upstream  temperature  T^  for  x  <  0  ,  and  an  arbitrary 
shear  stress  distribution  xQ(x,z)  for  x  _>.  0  .  The  amplitudes  of  the  step 
functions  e/h^f)  are  now  chosen  by  specifying 


(45) 


-11' 


1 


thus  equating  the  steps  to  increments  of  t^/h^  *  Doing  likewise  in  Equation 
(44)  and  superimposing  the  results,  we  obtain  the  inner  solution 


>2/3 


PC  a 
P 


1/3 


r  g(i»(  r 

£=0_  x1=5 


(46) 


for  the  case  of  an  adiabatic  boundary.  The  outer  limit  of  this  solution 
matches  the  inner  limit  of  Equation  (15).  The  outer  integrations  in  Eqs.  (45) 
and  (46)  are  begun  from  Z  =  0_  ,  with  tq  ini tally  zero  to  conform  to  the 
e -functions. 

Results  for  non-adiabatic  boundaries  may  be  obtained  by  combining  the 
present  results  with  those  of  Stewart  (1963).  In  that  work  Equation  (2)  was 
solved  without  the  dissipation  term,  with  a  uniform  upstream  temperature  Tm 
and  with  arbitrary  interfacial  thermal  conditions  for  x  >  0  .  Since  Eq.  (2) 
is  linear,  we  may  superimpose  Eq.  (25)  of  the  1963  paper  upon  the  particular 
solution  just  found;  this  gives  the  following  temperature  profile  for  the  case 
of  a  given  surface  temperature  distribution  TQ(x,z): 


T(x,y,z)  =  TA(x,y,z)  +J*  [1  -  n<n)]|^x  y  z  .  (47) 


Here  TQ  A  is  obtained  from  Eq.  (46)  with  y  =  0;  also 

‘l(n)  °  r(4'}3)  ^  exp^"  dn1  (48) 

with  n  defined  by  Equations  (24)  and  (32).  The  resulting  interfacial  heat 

flux  is 


.  3T| 

q0  3y|y=0 


k/hT 

z 


(49) 


/x  { r/v  wxii~1/3d(vTo,A)ic, 


1  /*»  J  l  J 

(9a) ,/Jr(4/3)  5=0_  x1=4 
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and  the  total  heat  transfer  rate  through  a  surface  region  0  <  x  <  x^(z)i 
Zj  <  z  <  zXI  is 

,1/3  z  x  (z)  x  (z) _ 

fi  -  - - - I  f  i  J  h  h  dx  )2/3d(T  -T  )|  .  (50) 

2r<4/3)a1/3  zx  5=0_  xf«t  Z  X  Z  1  0  °'A  t,Z 

A  two-dimensional  version  of  Eq.  (49)  was  given  by  Lighthill  (1950),  but  with 
an  approximate  function  TQ  based  on  the  calculations  of  Pohlhausen  (1921) 
for  Prandtl  numbers  near  unity.  The  present  results  are  preferable  for  large 
Prandtl  numbers. 

The  temperature  profile  corresponding  to  any  given  surface  heat  flux 
distribution  qQ(x,z)  is  obtained  by  combining  Eq.  (46)  above  with  Eq.  (39) 
of  Stewart  (1963): 

1  «  eXp(~  n3)q0hxhz  d4 

T(x,y,z)  »  T  ( x , y , z )  +  - - r - * - jjzz  )  . - - 57T  *  (51) 

1  r(l)  PCp(9o)Z/  0_  [/  -O  hxhzdXl)  ' 

This  completes  the  steady-state  solutions  to  lowest  order  in  a. 


A  -  .  4.-'  ' 


i 
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TRANSIENT  SOLUTIONS 


For  transient  processes,  an  outer  solution  for  a  ♦  0  is  obtained  by 

starting  from  Equation  (2)  without  the  conduction  term: 

-  DT(o) 

-  -  <T  :  VV)  .  (52) 

The  initial  condition  is  taken  to  be 


T(Xi,X2,X3,t)  -  T—  at 
Here  the  X^  are  material  coordinates, 
gives  the  time-dependent  outer  solution 
T(°)(X1,X2,X3,t)  -  T.  =  • 


t  *  0  .  (53) 

Integration  along  the  particle  paths 


PC  0 
P 


f  -(I  *  VV)| 


X  X  x  t  ($4) 


which  reduces  to  Equation  (15)  for  steady  flows. 

The  transient  inner  solution  for  the  dissipative  temperature  rise  is 
obtainable  from  Equation  (80)  of  Paper  II.  In  the  limit  of  small  a,  the 
dissipation  term  of  Eq.  (2)  may  be  approximated  in  the  boundary  layer  by  the 
local  interfacial  value.  The  resulting  inner  solution  for  T  is  independent 
of  the  normal  coordinate  y  ,  and  is  consequently  unaffected  by  heat 
conduction.  Equation  (54)  matches  this  adiabatic  inner  solution  exactly  at 
the  interface,  and  is  preferable  everywhere  else.  Thus,  for  adiabatic  time- 
dependent  systems  with  small  a  ,  Equation  (54)  is  valid  throughout  the 
fluid.  This  particular  integral  (with  T*°^  labelled  as  T^J  can  be 
superimposed  with  Eqs.  (76)-(79)  of  Paper  II  to  describe  unsteady-state  heat 
transfer  with  viscous  dissipation.  An  application  of  Eq.  (54)  is  given  in  the 
final  example. 
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APPLICATIONS  TO  SIMPLE  GEOMETRIES 


The  flows  considered  in  this  section  are  Newtonian  and  have  uniform 
temperature  upstream  from  a  locus  x  *  0  where  significant  dissipation 
begins.  Inner  solutions  are  presented  for  adiabatic  boundaries;  these  may  be 
superimposed  with  Equations  (49)  and  (51)  to  obtain  results  for  other 
interfacial  boundary  conditions. 


Developed  Flows  in  Ducts 

For  slow  flows  of  viscous  liquids  in  ducts,  the  hydrodynamic  development 
region  may  be  considerably  shorter  than  the  thermal  entrance  region.  The 
temperature  profiles  may  then  be  approximated  by  assuming  abrupt  flow 
development  at  the  entrance  of  the  duct.  Hie  x-coordinate  for  these  systems 
is  measured  parallel  to  the  flow  axis  and  the  z-coordlnate  is  measured  around 
the  duct  perimeter.  The  scale  factors  hx  and  hz  are  set  equal  to  unity. 
Several  duct  geometries  are  shown  in  Figure  3,  and  the  corresponding 
expressions  for  3  are  shown  in  Table  1  .  Equation  (46)  gives  the  inner 
solution  for  the  temperature  rise. 


(8iv2  Pr  x2  B4)1/3 


C 

P 


and  the  outer  solution  is  obtainable  from  Equation  (15). 


(55) 


Wedge  Flows 

Consider  a  laminar  two-dimensional  flow  impinging  on  the  wedge  shown  in 
Figure  4.  The  boundary-layer  coordinates  are  rectangular  with  hx  *  1  .  The 
velocity  profiles  are  given  by  Hartree  (1937);  the  interfacial  shear  rate  is 
3  =  f"(0,m)(  //2  (x  >  0)  (56) 
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and  a 


TABLE  1 


WALL  SHEAR  RATES  FOR  FULLY  DEVELOPED  NEWTONIAN 
FLOWS  IN  DUCTS 


Duct  Geometry 


Wall  shear  rate,  3 


Tube 


(P0  "  Pl)R/2u1, 


Annulus :  Inner  Cylinder 

Outer  Cylindet 

Equilateral  Three-walled  Duct 
(Landau  and  Lifshitz,  1960) 


IK 


Pl)r 


2  ML 


-K 


(1  +  <2) 

2  ln(l/ic) 


(Po  ~  pl>r  r  (,  -  <2> 

2ML  I  2  ln(1/K) 


2  ML 


Slit  with  One  Wall  Moving 
Downstream: 

Fixed  Wall 

Moving  Wall 


Pl)Y 


ML 


(Po- 


Pl)Y 


ML 


V 

w 

2Y 


Here  U  is  the  longitudinal  velocity  at  the  outer  edge  of  the  boundary  layer, 

U  =  ulXm  (  57 ) 

and  the  exponent  m  is  given  by  potential  flow  theory  as 

m  =  9/{ 2?r-8)  .  (58) 

The  dimensionless  wall  stress  fH(0,m)  is  tabulated  by  Hartree;  some  values 

are  f"(0,0)  *  0.4696  for  the  flat  plate  and  f"(0,1)  =  1.233  for  plane 

stagnation  flow.  Equation  (46)  gives  the  following  surface  temperature 

distribution  for  flow  over  an  adiabatic  wedge,  when  note  is  taken  of  the  jump 

in  8  (and  thus  in  t  /h  )  at  x  =  0  : 

0  2 


p(i) 

0,  A 


-  T  -  14 

p 


(59) 


For  the  special  case  of  a  flat  plate  (m  =  0),  the  resulting  recovery  factor 
is 


C  [t*1* 
Pl  0, A 


-  T-  J 


V2u 


»  1.92  Pr 


(60) 


in  agreement  with  the  asymptote  obtained  by  Meksyn  (196 0)  for  Pr  »  1  for 
this  geometry. 


Creeping  Plow  Around  a  Sphere 

Consider  the  steady  creeping  flow  of  a  Newtonian  fluid  past  a  stationary 
sphere.  The  surface  streamlines  are  contours  of  constant  longitude  ♦  , 
running  from  0  =  0  to  0  =  u  in  spherical  coordinates.  The  surface 
coordinates  for  the  inner  solution  are  x  *  0  and  z  *  $  ,  with  scale 
factors  hx  =  R  and  h^  =  R  sin0  .  The  Stokes  velocity  profile  gives  the 
interfacial  shear  rate 
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0 


(61) 


,  V  sin0 
3  oo 

2  R 

Thus,  t  /h  is  a  step  function,  vanishing  for  x  <  0  and  equal  to 
u  z 

2 

3mVj#/2R  for  x  >  0.  Equation  (46)  gives  the  dissipative  temperature  rise 

i  i  2 

T»  =  (  |  )3  9<<»  |  v-<  2  9  -  sin  29)  J3  (62) 

4  C  R 

P 

on  the  surface  of  an  adiabatic  sphere. 


(i 

0, 


Slow  Flow  Around  a  Cylinder 

Consider  the  slow  steady  flow  of  an  infinite  fluid  across  a  circular 
cylinder  of  radius  R  .  The  surface  streamlines  run  from  0=0  to  0  =  x  . 
The  surface  coordinates  for  the  inner  solution  are  x  =  9  and  z  =  Z  with 

t 

scale  factors  hx  =  R  and  hz  =  1,  where  (r,0,Z)  are  the  cylindrical 

coordinates.  A  series  expression  for  the  velocity  profile  near  the  cylinder 

is  given  by  Rosenhead  (1963).  Using  only  the  leading  terms,  one  obtains  the 

interfacial  shear  rate 

2V00  sin  0 

s  ■  — si -  <«> 

Here 

S  =  -  Y  +  In (8  Re-1 )  (64) 

and  Y  is  Euler's  constant.  Equation  (46)  gives  the  surface  temperature  rise 
due  to  viscous  dissipation 


2  4  1 

1296  Pr  v  v*  4 

(  - -  )3 

v  2  4  ' 

R  S 


g(0) 

T" 

c 

p 


1  2 

/**  (  J®  sin20  d9  )3  cos£d£ 

5=0_  ot-e  1 


(65) 


in  the  absence  of  heat  exchange  between  the  cylinder  and  fluid.  The  double 
integral  1(9)  has  been  evaluated  numerically  and  is  given  in  Table  2. 
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TABLE  2 


INTEGRAL  IN  EQ.  (65)  FOR  CYLINDER  SURFACE  TEMPERATURE 


180  e/n  1(0)  =  /9  (  /6  sin1/2e  de  )2/3cos&n 

£=o_  1 


* 

0 

0. 0000 

15 

0.0355 

30 

0.1377 

45 

0.2951 

60 

0.4896 

75 

0.6993 

90 

0.9008 

105 

1.073 

120 

1.198 

135 

1.267 

150 

1.278 

165 

1.243 

180 

1.195 

V 


^Stagnation  line. 


Flcv  Near  a  Rotating  Disk 


Consider  the  three-dimensional  axisymmetric  flow  induced  by  a  disk 
rotating  with  an  angular  velocity  u»  in  an  otherwise  quiescent  fluid.  A 
cylindrical  coordinate  system  is  chosen  which  rotates  with  the  disk.  From  the 
analysis  of  Cochran  (1934)  the  surface  velocity  components  are 

=  r  u  F (5)  (66) 

VQ  =  r  to  [G(U  -1]  (67) 

with 

C  =  Z  .  (68) 

From  Cochran's  tables  of  the  solution,  we  obtain 


Z=0 


2  2 

+  v2)2 


z=o 


(69) 


=  r  /o»3/V  /[F'  ( 0)  1  2  +  [G'  ( 0)  ]  2 

with  F'(0)  =  0.510  and  G'(0)  =  -0.616  . 

Since  the  temperature  is  independent  of  8  ,  we  choose  r  and  0  as 

surface  coordinates  (x,z)  for  the  inner  solution,  whence  hx  =  1  and  hz  = 

r.  This  is  simpler  than  integration  along  the  spiral  surface  streamlines,  and 

gives  the  same  result.  The  fundamental  solution  is  then  done  with  V0 

omitted,  and  with  0  replaced  by  0^  =  r  / a )3/v  F'(0)  [see  Equation  (69)] 

except  in  calculation  of  the  viscous  dissipation.  In  Equations  (45)  and  (46), 

we  must  accordingly  replace  e/h  0  by  e/h  0  ,  and  t  /h  by  T  0/h  8  • 

z  z  r  u  z  u  z  r 

The  resulting  inner  solution  for  the  dissipative  temperature  rise  is 


(wr)  (9  Pr ) 


1/3 


[ F' (0) ] 


2/3 


{ [F' (0) ] 2  +  [G‘ (0) ] 2}g 


3  /F'(O)a) 


3/2 


3  a  v 


1/2 


(70) 


for  a  thermally  insulated  disk. 
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Entrance  Flow  in  a  Tube 

Consider  a  fluid  which  enters  a  circular  tube  of  radius  R  with  an 
initially  flat  velocity  profile.  From  the  numerical  results  of  Hornbeck 
(1964),  the  quantity 


f  Re  = 


-4R 


3v 


<V>  3r 


r=R 


is  available  as  a  function  of  downstream  distance. 


(71) 


A  curve  fit  of  the  numerical  solutions  plotted  in  Figure  5  gives  the 
expression 


in  which 


f  Re 


1.328  q  2  +452  a1,1  +  15.6 

1  +  28.25  a1,1 


a  =  (x/R) (R<V>/v) 

Equation  (72)  satisfies  the  asymptotic  relations 


(72) 


(73) 


-  1 

lim  f  Re  «  1.328  a  2  (74) 

a+0 


lim  f  Re  -  16  .  (75) 

ax* 

Equation  (46),  (71)  and  (72)  give  the  surface  temperature  rise  in  the  thermal 
entrance  region  of  an  adiabatic  tube: 


0,  A 


T  = 


(tI  )3g(°) 


Pr3  <V>2 


r 

C=o 


(  J  /  f  Re  d  a  ) 3  d(f  Re) | 

_  r  » 


V* 


(76) 

This  solution  is  plotted  in  dimensionless  form  in  Figure  6.  Near  the 
entrance,  the  temperature  rise  is  essentially  that  of  the  flat  plate,  whereas 
downstream  the  curve  approaches  the  asymptote  for  hydrodynamically  developed 
laminar  flow. 
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Wall  shear  stress  function  for  entrance  flow  in  a  tube. 


1 


A  WIRE-COATING  PROBLEM 

Improvements  in  wire-coating  processes  have  been  made  largely  by 
empirical  means.  Experimental  studies  of  such  flows  are  difficult,  because  of 
the  small  clearances  required  in  a  wire-coating  die.  Theoretical  treatment  is 
hampered  by  temperature-dependent  properties,  along  with  complicated  rheology 
and  geometry.  Previous  theoretical  work,  summarized  by  Carley,  End<b,  and 
Krantz  (1978),  has  included  approximate  analytical  treatments  and  elaborate 
numerical  computations.  In  this  section  we  give  an  asymptotic  analysis  of  a 
representative  wire-coating  process. 

The  geometry  (see  Figure  7)  and  operating  conditions  of  a  typical  wire¬ 
coating  die  are  given  by  Haas  and  Skewis  (1974).  The  die  wall  is  tapered, 
with  local  radius 

R2(x)  =  a(x)R1.  (77) 

Here  R1  is  the  wire  diameter,  and 

a(x)  -  a(0)  (1  -  ~  )  +  a (L) (  £  )  (78) 

with 

a(0)  «  3.39 
a(L)  *  1.17 

L  »  0.39  in.  (0.99  cm) 

R,  *  0.0127  in.  (0.0323  cm) 

The  coordinate  x  is  the  axial  distance  from  the  inlet  of  the  die.  The  wire 
velocity  vw  is  4000  ft.  min-1  =  2032  cm  S"1. 

An  exact  description  of  the  flow  field  would  require  numerical  integra¬ 
tion  of  the  coupled  mass,  momentum,  and  energy  equations.  Here  an  approximate 
velocity  profile  is  used.  The  axial  velocity  profile  in  the  melt  is 

approximated  by  a  three-term  expansion 
V 

-2-  =  f  (x)  +  f ,  (x )ln  (J-  )  ♦  f  <x)(jj-  V  (79) 

w  1  1 
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based  on  functions  appropriate  to  annular  Newtonian  flow  with  a  moving  inner 
cylinder.  Use  of  the  no-slip  conditions  at  the  interfaces# 


V  =  V 
x  w 


at  r  =  R, 


Vx  =  0  at  r  =  a(x)R1 


(80) 
(81  ) 


and  the  given  mass  flow  rate, 


aR. 


1 

w  =  2np  /  V  rdr 


(82) 


leads  to  three  equations  in  the  three  unknowns: 

f0(x)  +  f2(x)  =  1 


f q ( x)  +  fj(x)  In  a(x)  +  f2(x)a  (x) 


(83) 

I  (84) 

■  A. 

a*(L)  -  1  =  /aU)  2 (fQ ( x)  +  f^x)  In  u  +  f2(x)u2)udu  .  (85) 

Here  w  has  been  chosen  so  that  there  is  no  draw-down,  i.e.,  the  final 
coating  diameter  matches  the  exit  diameter  of  the  die. 

From  Equation  (79)  the  shear  rates  at  the  die  wall  and  the  wire  surface 


are 


0d  = 


V  f  (X) 

r  (  ZuT  *  2‘,<’0£2(’°] 


(86) 


0w  = 


(87) 


57 

These  functions  are  plotted  in  Figure  8.  The  shear  rate  at  the  die  wall  rises 


monotonically  to  a  maximum  of  3,3  *  105  sec  1 


wire  first  declines  to  zero  and  then  rises  to  2.1  x  IQ5  sec  1 


while  the  shear  rate  at  the 
In  the 

entrance  region,  the  fluid  velocities  do  not  exceed  the  wire  velocity.  In  the 
exit  region,  however,  fluid  velocities  greater  than  the  wire  velocity  must 
occur  to  satisfy  the  zero  draw- down  specification. 

Carley,  Endo,  and  Krantz  (1978)  have  provided  property  data  for  a 
representative  low-density  polyethylene  melt.  For  the  boundary  regions,  they 
use  the  shear  stress  model 
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"  4  V 
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T 
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dV 
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dr 


ref 


mQ  =  101,000 

E  *  11,700 
R  »  1.987 


g  cm 


-1 


sec 


-5/: 


cal  mole 
cal  mole  1  K 


ref 


479  K  . 


(88) 


(89) 


In  this  analysis  m^  is  taken  as  a  constant,  and  is  evaluated  as  described 
below.  The  following  properties  are  treated  as  independent  of  temperature: 

k  =  0.000786  cal  sec  1  cm  1  K  1 
p  =  0.91  g  cm-3 
Cp  *  0.55  cal  g-1  K-1  . 

The  surfaces  of  the  wire-coating  system  are  assumed  adiabatic.  At  the 
die  wall  we  use  hx  =  1  and  h2  =  a(x)R1  .  Equations  (46),  (86),  and  (88)  then 
give  the  die  wall  temperature  distribution 


TU)  -  T 

0  ■ 


9  2/3  g(  0) 


K  0;  (  rx 

1/3  mT  '  rt 
o'  x^*0 
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*  f  (f 
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for  adiabatic  operation.  The  first  term  of  Equation  (90)  accounts  for  the 
assumed  initial  jump  in  wall  shear  rate,  and  the  second  term  accounts  for  all 
subsequent  changes.  The  index  mT  is  evaluated  from  Equation  (89)  at  the 
mean  of  the  entrance  and  exit  surface  temperatures;  this  entails  an  iterative 
calculation. 

The  temperature  rise  at  the  moving  wire  surface  is  obtained  from 
Equations  (54),  (87)  and  (88),  with  t  *  x/Vw  and  -(T  :  Vv)  *  TQ3w.  The 

result  is 


T^°^  —  T 
0  « 


PC  V 
p  w 


r .  ».3 «,  • 


The  index  m^  is  again  evaluated  at  the  mean  of  the  entrance  and  exit  surface 
temperatures.  Equation  (15)  is  equivalent  to  Eq.  (54)  here  since  the  system 
is  at  steady  state  when  viewed  from  the  die. 

Numerical  integration  of  Equations  (90)  and  (91)  with  T—  *  Tref  9*ves 
the  results  shown  in  Figure  9  .  The  temperature  at  the  die  wall  rises  from 
479K  to  534K;  this  increase  is  comparable  to  the  results  obtained 
numerically  by  Car ley,  Endo,  and  Krantz  (1978)  for  similar  wire-coating 
systems.  The  calculated  temperature  rise  at  the  wire  surface  is  much  smaller 
(only  5.8K  at  the  maximum)  and  becomes  still  less  if  one  allows  for  heat 
transfer  to  the  moving  wire. 


Figure  9.  Surface  temperature  rises  for  an  adiabatic  wire-coating  operation 


The  boundary  layer  thickness  associated  with  each  increment  of  T0/^z 


; 

alonq  the  die  wall  is  obtainable  from  Equations  (32)  and  (86): 

<5(x,5)  = — - -  {9a  f*  /~afl7  a  ^x.}3  (92) 

J  rt  f  U  1 

Xi=,  9 

Piqure  (10)  shows  the  qrowth  of  five  representative  boundary  layers  emanatinq 
from  equally  spaced  points  alonq  the  die  wall.  Each  of  these  layers  qrows 
rapidly  at  first,  but  then  qets  thinner  as  the  flow  converqes  and  acceler¬ 
ates.  The  final  thickness  of  the  layer  associated  with  the  initial  jump  in 
tq  is  0.0011  inches  (0.0028  cm)  or  11%  of  the  final  coatinq  thickness. 

Thus,  the  small-diffusivity  assumption  of  the  analysis  is  accurate  here. 
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NOTATION 

=  lenqth  of  side  of  equilateral  triangle?  radius  ratio  in  Eqs. 

( 77)-( 78). 

=  P(m)r (n)/r(m+n) 

=  heat  capacity  at  constant  pressure  (EM-1T-1) 

)  > 

=  ~Jr  +  (V  *  V),  substantial  derivative 
=  dimensionless  function  defined  in  Equation  (66) 

=  friction  factor,  see  Equation  (71) 

=  function  defined  in  Equation  (23)  (T) 

=  dimensionless  shear  stress  at  wedge  surface,  see  Equation  (56) 
=  dimensionless  functions  defined  in  Equation  (79) 

=  dimensionless  function  defined  in  Equation  (67) 

=  dimensionless  temperature  function  defined  in  Equation  (23) 

=  scale  factors  defined  in  Equation  (3) 

*  thermal  conductivity  (Et-1^1!-1 ) 

=  lenqth  of  channel 
=  arc  length  alonq  a  streamline 

=  dimensionless  parameter  defined  in  Equation  (58) 

=  parameter  in  power  law  model,  see  Equation  (88) 

=  terms  of  order  n  or  qreater  in  y 
=  p  +  pQ  (ML-1t~2) 

=  decrease  of  P  alonq  duct  length  L 
=  Prandtl  number 
=  pressure  (M  L-1t~2) 

=  conductive  heat  flux  into  the  stream  at  y  =■  0  (E  L-2t-1 ) 


•  *■!* 
W  ' 


R  =  gas  constant  (E  mole-1T-1 ) 

R  =  radius  of  sphere  or  cylinder 

Re  =  Reynolds  number 

R1  =  wire  radius 

r2(x)  =  local  radius  of  wire-coating  die 

r  =  position  vector 

r  =  radial  coordinate 

r  =  dimensionless  ratio  C^/Cq  »  see  Equation  (34) 

T  =  absolute  temperature 

U  =  local  free  stream  velocity  for  wedge  flow  (Lt-1 ) 

<V>  =  mean  velocity  through  cross-section  of  pipe  (Lt-1) 

V  =  velocity  vector  relative  to  non-deforming  coordinates  (Lt-1 ) 

V  =  magnitude  of  velocity  vector  (Lt-1 ) 

vw  =  velocity  of  wire  or  moving  wall  (Lt-1 ) 

Vx,Vy,Vz  =  components  of  velocity  vector  (Lt-1 ) 

V  ,V„  =  velocity  components  in  cylindrical  or  spherical  coordinates 

r  o 

(Lt-1) 

=  tangential  component  of  velocity  vector  (Lt-1 ) 

=  magnitude  of  free  stream  velocity  (Lt-1 ) 

X1fX2,X3  =  material  coordinates 

x  =  surface  coordinate  directed  along  surface  streamlines; 

see  Equations  <3)-(6)  (Lhx-1) 

V  =  half-width  of  slit  (L) 

y  =  normal  distance  from  interface  (L) 

Z  =  cylindrical  axial  coordinate  (L) 

z  =  surface  coordinate  normal  to  the  local  x-direction; 

see  Equations  <3)-(6)  (Lh  “1 ) 
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Greek  Symbols 

a  *  k/pC  ,  thermal  diffusivity  (L2t_1) 

P 

8  =  (3vt/3y)|  g,  function  introduced  in  Equation  (4)  ( t—  1 ) 

r (x)  =  J  e  2  zX  ^dz 

0 

Y  =  Euler's  constant,  0.57721... 

6  *  boundary  layer  thickness  for  viscous  heating  or  heat  transfer 

<L) 

6  ,5  ,6  =  unit  vectors 

~x  ~y  ~z 

e  =  special  dissipation  function  for  fundamental  solution,  see 

Equations  (19)  and  (43) 

n  3  dimensionless  coordinate  defined  in  Equation  (24) 

0  =  angle  in  cylindrical  or  spherical  coordinates 

0  =  wedge  angle  (radians) 

k  =  ratio  of  inner  cylinder  radius  to  outer 

M  *  viscosity  (ML~^t_1 ) 

v  *  p/p  kinematic  viscosity  (L2t-1 ) 

S  *  value  of  x  at  the  jump  in  e  and 

K  =  dimensionless  temperature  function  in  Equation  (48) 

P  =  density  (ML~^) 

o  =  (x/R) (R  <V>/v)  1  , dimensionless  downstream  coordinate  in  tube 

t  =  stress  tensor  (ML_1t*2) 

tq  *  magnitude  of  tangential  stress  at  interface  (ML_1t~2) 

A 

♦  *  gravitational  potential  energy  (L2t*2) 

$  *  angle  in  spherical  coordinates 

w  =  angular  velocity  of  rotating  disk  (radians  sec*1) 

♦  *  stream  functions,  see  Equation  (7) 
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Superscripts 


(i)  inner  solution 

(o)  outer  solution 

*  solution  with  special  dissipation  function 


Subscripts 


A  dissipative  solution  for  adiabatic  boundary 

d  at  die  wall 

w  at  wire  surface 

0  in  or  into  the  fluid  at  y  =  0 

•  upstream  or  initial  state 

1  dummy  variable 


e(x,£,z 
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